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j> . ABSTRACT 

QO ' We present an exhaustive analysis of five broad-band observations of GRS 1915+105 

t— I , in two variability states, x an d u>, observed simultaneously by the PCA and HEXTE 

■ detectors aboard RXTE and the OSSE detector aboard CGRO. We find all the spec- 
tra well fitted by Comptonization of disc blackbody photons, with very strong ev- 
idence for the presence of a nonthermal electron component in the Comptonizing 

■ plasma. Both the energy and the power spectra in the x state are typical to the very 

high/intermediate state of black-hole binaries. The spectrum of the lj state is char- 

<i ■ acterized by a strong blackbody component Comptonized by thermal electrons and a 

Oh' weak nonthermal tail. We then calculate rms spectra (fractional variability as func- 

q \ tions of energy) for the PCA data. We accurately model the rms spectra by coherent 

5^ ■ superposition of variability in the components implied by the spectral fits, namely a 

' less variable blackbody and more variable Comptonization. The latter dominates at 

■ high energies, resulting in a flattening of the rms at high energies in most of the data. 
^ ■ This is also the case for the spectra of the QPOs present in the x state. Then, some of 

our data require a radial dependence of the rms of the disc blackbody. We also study 
the distance to the source, and find d ~ 11 kpc as the most likely value, contrary to 
' a recent claim of a much lower value. 
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1 INTRODUCTION 

GRS 1915+105 is a binary system consisting of an early- 
type K giant (Greiner et al. 2001b) accreting via a Roche- 
lobe overflow onto a black hole with the mass function of 
9.5±3M Q (Greiner, Cuby & McCaughrean 2001a). The re- 
cent measurement of the rotational broadening of the giant 
combined with the orbital parameters has resulted in the de- 
termination of the black-hole and donor mass of 14 ± 4Mq 
and 0.8 ± 0.5M Q , respectively (Harlaftis & Greiner 2004). 
The accretion onto the black hole yields highly variable 
X-ray emission (e.g., Belloni et al. 2000, and references 
therein). Still, even the hardest observed spectra are rela- 
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tively soft, consisting of a blackbody-like component and a 
high-energy tail (e.g., Vilhu et al. 2001). They are softer than 
those of other black-hole binaries in the hard state, which 
EFe spectra peak at ~ 100 keV (e.g., Cyg X-l, Gierlihski 
et al. 1997), but generally similar to their intermediate and 
soft states (e.g., XTE J1550-564, Gierlihski & Done 2003; 
Cyg X-l, Gierlihski et al. 1999; LMC X-l, LMC X-3, Wilms 
et al. 2001). This is confirmed by the comprehensive analy- 
sis of the states of GRS 1915+105 by Done, Wardzihski & 
Gierlihski (2004), who attribute the lack of the hard state 
in GRS 1915+105 to its very high, Eddington or higher, 
luminosity. 

The blackbody component arises, most likely, in an 
optically-thick accretion disc. On the other hand, there were 
some uncertainty regarding the origin of the tail. Three main 
models proposed involved Comptonization of the blackbody 
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photons by high-energy electrons. They differed, however, 
in the distribution (and location) of the electrons, which 
are assumed to be either thermal (Maxwellian) , non-thermal 
(close to a power law), or in a free fall onto the black hole. 

A discussion of these models is given in Zdziarski (2000) , 
who shows that the thermal and free-fall models of the soft 
state of black hole binaries can be ruled out, mostly by the 
marked absence of a high-energy cutoff around 100 keV in 
the CGRO data (Grove et al. 1998; Gierlihski et al. 1999; 
Tomsick et al. 1999; Zdziarski et al. 2001, hereafter Paper 
I; McConnell et al. 2002). In particular, those two models 
have been shown to be strongly ruled out for GRS 1915+105 
in Paper I. The present best soft-state model involves elec- 
tron acceleration out of a Maxwellian distribution (i.e., a 
non-thermal process), which leads to a hybrid electron dis- 
tribution consisting of both thermal and non-thermal parts 
(Zdziarski, Lightman & Maciolek-Niedzwiecki 1993; Pouta- 
nen & Coppi 1998; Gierliriski et al. 1999; Coppi 1999). 

In order to gain further insight into the nature of ra- 
diative processes in GRS 1915+105, we study here avail- 
able broad-band, ~2-10 3 keV, spectra and variability of 
this source from simultaneous observations by the PCA and 
HEXTE detectors on board RXTE and the OSSE detec- 
tor on board CGRO. An earlier study of two of those spec- 
tra corresponding to the lowest and highest X-ray flux was 
given in Paper I. The spectra, showing extended power laws 
without any cutoff up to at least 600 keV, provide strong 
evidence for the presence of non-thermal Comptonization. 
Here, we also study the rms spectra from the PCA data and 
their implications for physics of the source. Given the im- 
portance of the distance and inclination of the system for 
our spectral fits and their theoretical interpretation, we be- 
gin our analysis with estimating those parameters in Section 
H below. 



2 THE DISTANCE AND INCLINATION 

Using the rotational velocity of the donor, v sin i — 26 + 3 km 
s _1 (Harlaftis & Greiner 2004), and assuming synchronous 
rotation with the orbital period, P or b = 33.5 d, we find the 
donor radius of R2 = (17.2 + 2.0) sin -1 i Rq. The reddening- 
corrected secondary magnitude is Kq = 12.1-12.6, as ob- 
tained using the unveiled magnitude of the donor, K — 14.5- 
15 (Greiner et al. 2001b), and correcting it for the extinction 
of A K = 2.4 ± 0.2 (Chapuis & Corbel 2004). Then, we use 
the Barnes- Evans relation (Cahn 1980), 

Fk = 4.2211 — 0.1K — 0.5 log 10 s, (1) 

where the K surface brightness is Fk = 3.84 ( Beuermann, 
Baraffe & Hauschildt 1999), which yields the angular di- 
ameter of the secondary of s = 0.0220-0.0175 mas. This 
corresponds to the distance of 

d = (6.4-10.2) sin _1 i kpc, (2) 

or d = 6.9-11.3 kpc for i = 66 ± 2° (Fender et al. 1999; see 
also Mirabel & Rodriguez 1994). 

Other estimates include 6 < d < 12 kpc from the kine- 
matics of interstellar clouds (Dhawan Goss & Rodriguez 
2000a; Chapuis & Corbel 2004), d < 11.2 ± 0.8 kpc from 
radio data (Fender et al. 1999), d ~ 12 kpc from the secular 
(transverse) motion (Dhawan, Mirabel & Rodriguez 2000b), 



and, the most accurate, distance of d = 12.1 ± 0.8 kpc im- 
plied by the radial systemic velocity (Greiner et al. 2001a). 
The mutual agreement of the last two estimates is consis- 
tent with the motion of the source being primarily due to 
the Galactic rotation. 

Recently, Kaiser et al. (2004) have identified two IRAS 
sources, IRAS 19124+1106 and IRAS 19132+1035, each lo- 
cated at the angular distance of ~ 26' on opposite sides of 
the direction to GRS 1915+105, with the impact sites of 
the jets of that source. This identification resulted in a new 
distance estimate to GRS 1915+105, d ~ 6.5 kpc, which, 
in combination the observed jet motions on small scales, 
yielded a jet velocity of about 0.7c and an angle of 47-53° 
of the jets to our line of sight. However, the two IRAS re- 
gions do not move together with GRS 1915+105. The ra- 
dial velocities of the two IRAS sources are redshifted by 
about 70 km s _1 (IRAS 19132+1035) and 60 km s _1 (IRAS 
19124+1106), respectively, with respect to the radial sys- 
temic velocity of GRS 1915+105 of -3+10 km s _1 derived 
from studies of 12 CO and 13 CO lines in its K-band spec- 
tra (Greiner et al. 2001a). Kaiser et al. (2004) doubted the 
absolute wavelength calibration of the observed spectra and 
suggested that 'systematic effects may well have altered the 
determination of the systemic velocity'; however they nei- 
ther clarified nor proposed any such effects. On the other 
hand, it is known from studies of symbiotic red giants that 
the radial velocity curves based on near-IR CO lines are very 
accurate, and any systematic effects are < 1-2 km s~ (e.g., 
Fekel et al. 2000). In particular, there is an excellent agree- 
ment between the systemic velocity derived from CO and 
optical lines. Although in the case of GRS 1915+105, the 
CO band accuracy is less than in typical sources due to rel- 
atively high veiling of the giant spectrum, it is unlikely that 
the error exceeds about 8 km s _1 (J. Greiner, personal com- 
munication). Moreover, this lower distance is incompatible 
with the lower limit derived from the red giant radius and 
brightness, which for the system inclination (= jet inclina- 
tion) of i = 53° yields, using equation @, d — 8.0-12.8 kpc, 
i.e., substantially more than the distance used by Kaiser et 
al. (2004) to obtain that inclination. 

Thus, the two IRAS sources appear not to be associated 
with GRS 1915+105. Combining other constraints above, we 
find d ~ 11 kpc as the only value approximately compatible 
with all of them. Hereafter we assume this distance and i — 
66°. 



3 THE DATA 

We use RXTE data simultaneous with those by the OSSE 
and showing relatively weak X-ray variability, which allows 
a meaningful analysis of the average spectra. In Paper I, the 
variability states used for spectral fitting were called x an d 
7, according to the classification of Belloni et al. (2000). 
However, Klein- Wolt et al. (2002) introduced then an ad- 
ditional variability state, to, and Naik, Rao & Chakrabarti 
(2002) identified the 7-state observation of Paper I as be- 
longing to the state uj, which name we use hereafter. 

We extract spectra from pointed RXTE observations 
using FTOOLS v. 5.3. We add systematic errors of 1 per cent 
to each PCA and HEXTE channel. We extract spectra from 
both HEXTE clusters separately. Since those data sets were 
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obtained with various number of PCA detector units (PCUs) 
operating, we have decided to use only the PCU and 2, 
which were switched on in all of our data. We use only layer 
1 of the PCA. In order to test the effect of changing the 
method of data extraction, we have also obtained the PCA 
spectra for all operating detectors for a given data set, and 
used all layers. However, we have found only negligible dif- 
ferences in the results of spectral fits to those data with 
respect to the PCU 0, 2, layer 1, configuration. We have 
also tested for the effect of the choice of a particular RXTE 
observation in cases when more than one exist (for the \ 
state, see Paper I) , which effect has been found to be minor. 
Table [U shows the log of the RXTE and OSSE data used. 

Since we need to integrate over rather long periods 
in order to obtain accurate spectra, we have searched for 
flares/dips from the lightcurves prior to the spectral extrac- 
tion. This affected the x-state RXTE spectrum of VP 619, 
were we removed a few flare-like events. Furthermore, the 
lightcurve of the ui state, VP 813, consists of steady emis- 
sion regularly interrupted by short dips. For that spectrum, 
we have then obtained separately the RXTE spectra for the 
steady state and the dips, as shown in Fig. We note that 
the steady emission outside the dips is very similar to that 
of the class 7, which classification was used in Paper I. 

The OSSE detector accumulated spectra in a sequence 
of 2-min measurements of the source field alternated with 2- 
min, offset-pointed measurements of background. The back- 
ground spectrum for each source field was derived bin-by- 
bin with a quadratic interpolation in time of the nearest 
background fields (Johnson et al. 1993). The OSSE data 
include energy-dependent systematic errors estimated from 
the uncertainties in the low-energy calibration and response 
of the detectors using both in-orbit and pre-launch cali- 
bration data. They range from 0.03 at 50 keV to 0.003 at 
> 150 keV (where, however, statistical errors dominate). The 
OSSE data used here are the same as in Paper I. 

Apart from the energy spectra, we also calculate the 
corresponding power density spectra (PDS) from PCA 
lightcurves with 1/128-s resolution. Each PDS is an average 
of several Fourier transforms calculated over 512-s intervals 
of the lightcurve. This gives the frequency band of 1/512-64 
Hz. We create PDS both in the full 2-60 keV PCA energy 
band, and in three bands of approximately 2-7, 7-15 and 
15-60 keV. 

In order to study energy dependence of variability in de- 
tails, we also extract energy-dependent rms, using two differ- 
ent methods. In both methods, we first extract background- 
subtracted lightcurves with 1/128-s resolution in each PCA 
energy channel. In the first method, we calculate intrinsic or 
excess rms [see eq. (10) in Vaughan et al. 2003], characteriz- 
ing the intrinsic source variability not related to the Poisso- 
nian noise. In the second method, we calculate the PDS from 
each of the lightcurves (over 512-s intervals), subtract the 
Poissonian noise, and correct for dead-time effects (Revnivt- 
sev, Gilfanov & Churazov 2000). The energy-dependent rms 
is then found by integrating the PDS over 1/512-64 Hz fre- 
quency band. Both methods give very similar results, see 
Section below. 
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Figure 1. The PCA lightcurve of the observation in the u state 
(VP 813), showing the intervals used for calculating the high and 
low spectra. The origin of the time axis corresponds to 1999-04-23 
02:05:03 UT. 

4 SPECTRAL MODELS 

We first fit the data with the hybrid Comptonization model, 
COMPPS, of Poutanen & Svensson (1996), as implemented in 
XSPEC (Arnaud 1996). In this model, the temperature, kT, 
and the total Thomson optical depth, r, of the electrons 
are free parameters and the distribution is assumed to be 
Maxwellian up to the Lorentz factor, 7 m in (which is the third 
free parameter), and a power law with an index, Y, above it 
up to 7max, assumed here to be 10 3 . The iterative scattering 
method of Poutanen & Svensson (1996) is not applicable 
to plasmas with large t, and we constrain t to < 6. We 
assume a spherical geometry of the hot plasma (the model 
geometry parameter is set to 0). Since some of the data may 
only weakly constrain the form of the nonthermal tail, we 
assume F > 2, which does not allow for very hard tails, 
unlikely to be produced by a physical acceleration process. 
The seed photons for Comptonization are assumed to have 
the spectrum of a blackbody disc (Mitsuda et al. 1984) with 
the maximum temperature of kTbb- 

We also take into account Compton reflection with a 
solid angle of Q (Magdziarz & Zdziarski 1995) and an Fe 
Kq emission from an accretion disc (Shakura & Sunyaev 
1973) assumed to extend down to 100GM/c 2 . We allow the 
reflecting surface to be ionized, with the ionization param- 
eter of £ = Lion/nR 2 , where L- lon is defined as the 5 eV-20 
keV luminosity in a power law spectrum and n is the den- 
sity of the reflector located at distance R from the illumi- 
nating source. We use here the ionization model by Done 
et al. (1992), which is not applicable to very highly ionized 
Fe (Ballantyne, Ross & Fabian 2001), but sufficiently ac- 
curate at low/moderate ionization, such as that found in 
our data. Consistent with the found low ionization, we keep 
the line energy fixed at 6.4 keV. The reflector temperature 
is kept at 10 6 K. The elemental abundances of Anders & 
Ebihara (1982), an absorbing column of Ah > N§ , and 
i — 66° (Section |5J| are assumed. Here, N§ is the Galactic 
column density in the source direction, which we take here as 
~ 1.8x 10 22 cm" 2 , following Dickey & Lockman (1990). (We 
note, however, that a higher value of N§ has been claimed 
by Chaty et al. 1996.) 

We then fit the data with the model eqpair 
(Coppi 1999; Gierliriski et al. 1999), which calculates self- 
consistently microscopic processes in a hot plasma with elec- 
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Table X. Log of the OSSE and RXTE observations. 





OSSE 






RXTE 




VP 


Start 


End 


Observation ID 


Date 


Time (UT) 


Variability state 


601 


1996-10-15 


1996-10-29 


10408-01-45-00 


1996-10-29 


11:56-17:30 


X 


619 


1997-05-14 


1997-05-20 


20187-02-02-00 


1997-05-15 


11:31-18:42 


X 


720 


1998-05-05 


1998-05-15 


30402-01-12-01 


1998-05-12 


00:28-01:00 


X 


813 


1999-04-21 


1999-04-27 


40403-01-07-00 


1999-04-23 


01:57-02:44 


UJ 


917 


2000-04-18 


2000-04-25 


50405-01-02-00 


2000-04-22 


09:21-10:19 


X 



tron acceleration at a power law rate with an index, Fi n j, in a 
background thermal plasma with a Thomson optical depth 
of ionization electrons, n. The electron temperature, kT, 
is calculated from the balance of Compton and Coulomb 
energy exchange, as well as e pair production (yielding 
the total optical depth of r > n) is taken into account. 
The last two processes depend on the plasma compactness, 
I = Cot I \lZm c c 3 ) , where C is a power supplied to the hot 
plasma, 1Z is its characteristic size, and or is the Thomson 
cross section. We then define a hard compactness, 4», corre- 
sponding to the power supplied to the electrons, and a soft 
compactness, £ s , corresponding to the power in seed photons 
(of a disc blackbody spectrum) irradiating the plasma. We 
assume £ s — 100, which is consistent with the source lumi- 
nosity and estimates on the size of the hot plasma, see Paper 
I. The compactnesses corresponding to the electron acceler- 
ation and to a direct heating (i.e., in addition to Coulomb 
energy exchange with non-thermal e and Compton heat- 
ing) of the thermal e are denoted as £ nt h and £ t h, respec- 
tively, and £h — £ n th + £th- Details of the model are given in 
Gierlihski et al. (1999). 

The nonthermal electrons are accelerated between the 
Lorentz factors 7 m i n and 7 ma x, with the latter assumed to be 
10 3 (as in the compps model). For 7 m i n , we use the best-fit 
value of the COMPPS, as 7min in eqpair only weakly affects 
the fit results and thus cannot be reliably fitted. Reflection 
and absorption is treated the same way as in COMPPS. 

We note that results from the CG720/EGRET, which 
did not detect GRS 1915+105, yield an upper limit at ~100 
MeV of EF E <0.04: keV cm" 2 s" 1 (Hartman et al. 1999). 
This upper limit is fully satisfied in all of eqpair fits to the 
data obtained by us below, since that model takes into ac- 
count photon absorption in pair-producing photon-photon 
collisions. Also, that model treats self-consistently energetic 
electrons and positrons produced in such events. On the 
other hand, the COMPPS model does not include pair pro- 
duction, and some of our fit results may violate the EGRET 
upper limit. Thus, our results obtained with COMPPS should 
be treated as applicable only to the photon energy range up 
to soft 7-rays, and our preferred overall model is eqpair. On 
the other hand, COMPPS treats Comptonization more accu- 
rately as well as it gives us directly the steady-state electron 
distribution in the source. 

The uncertainties ranges for each parameter below are 
given for 90 per cent confidence, i.e., A\ 2 = 2.71. To as- 
sess the significance of the nonthermal electrons, we also 
fit the data with purely thermal models (i.e., without any 
high-energy tail), using both compps and eqpair. We then 
use the F-test (Bevington & Robinson 1992) to calculate the 
probability, P t h, that the fit improvement by allowing a non- 




E [keV] 

Figure 2. Broad-band spectra in the \ state. The solid and 
dashed curves give the intrinsic and absorbed EQPAIR model spec- 
tra, respectively. These and subsequent spectra are shown normal- 
ized to the PCA data. 

thermal electron tail (with two additional free parameters, 
the index and normalization) is by chance. 



5 BROAD-BAND SPECTRA AND PCA 
POWER SPECTRA 

The fit results for both COMPPS and eqpair models are given 
in Table 2. We find rather good agreement between those 
results, in particular for the blackbody and electron tem- 
peratures and the absorbing column and the characteristics 
of reflection. The values of r are somewhat different as they 
are defined in different way in the two models. Note that in 
most cases F > r^j , due to the softening of the steady-state 
electron spectrum (fitted in compps) with respect to the 
spectrum of the injected electrons (fitted in eqpair). Also, 
the results for VP 619 and 813 are similar to those given in 
Paper I. In all cases, the values of the relative normaliza- 
tion between the PCA, HEXTE and OSSE spectra are very 
close to unity. The values of \v < 1 are due to the systematic 
errors being apparently somewhat overestimated. 

The bolometric isotropic fluxes given in Table 2 corre- 
spond to L/L E » 7x lO 6 F b oi(-D/llkpc) 2 (14M /M), where 
L E « 1.5(M/Mq) x 10 38 erg s" 1 is the Eddington luminos- 
ity. They yield L/L E « 0.17-0.24 and 0.5-0.6 for the x states 
and the high-flux part of the u state, respectively. 

We confirm the conclusion of Paper I that the 
commonly-used model of a disc blackbody and an e-folded 
power law provide very bad fits to all of our data, with 
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Figure 3. Spectral components of the EQPAIR fits (solid curves) 
to the x state. The dots show the unscattered blackbody, the 
short and long dashes show Compton scattering by thermal and 
nonthermal electrons, respectively (the latter also includes the 
pair annihilation component), and the dot-dashed curves show 
Compton reflection including the Fe Ka line. 



X v 2> 1- On the other hand, all our fits with the physi- 
cal Comptonization model yield xt < 1) see Table 2. All 
the data rule out purely thermal Comptonization, with ex- 
tremely low values of the corresponding probability, P t h- Be- 
low, we discuss our results separately for the x an d u states. 




Figure 4. The distribution of the Comptonizing electrons cor- 
responding to the best fit with the COMPPS model to the x-state 
data of the VP 619. See Table 2 for the parameters. The dis- 
tribution of the momentum, TV(lnp) = pN(p) is shown, and its 
normalization is arbitrary. The distibution is Maxwcllian up to 
the p corresponding to the 7 m i n and a power law, <x 7~ r above 
it. 



5.1 The x state 

Fig. |21 shows the unfolded spectra and the best-fit eqpair 
models for our four data sets in the x state. We see all the 
spectra are rather similar, with relatively modest variability 
on timescales between the observations. Two pairs of spectra 
differ by the level of the ~10 keV flux (high: VP 601, 720, 
low: VP 619, 917) but there is no obvious correlation with 
the level of the high-energy, ~100 keV flux. 

The two spectra with the high 10-keV flux also show 
relatively strong reflection, see Table 2. The spectral com- 
ponents for the four spectra are shown in Fig.|3] Fig.^Jshows 
the electron distribution corresponding to the COMPPS best- 
fit model to the spectrum of the VP 619. 

The thermal Compton component shown in Fig.[3]corre- 
sponds to emission from purely thermal plasma of the same 
r and kT (adjusted in eqpair by a suitable reduction of 
£h/£s) as the thermalized electrons in the hybrid distribu- 
tion. The non-thermal component is the difference between 
the total scattering spectrum by the hybrid electrons and the 
thermal-scattering spectrum. See Hannikainen et al. (2005) 
for details of the method. 

We see in Fig. [3] that our models predict the power- 
law emission extending with no cutoff up to several MeV 
and weak annihilation feature (with the plasma being pair- 
dominated in some cases, e.g., with n — » allowed for VP 
619 and 917). These predictions should be tested by future 
soft 7-ray missions. On the other hand, we predict the emis- 
sion at ~100 MeV to be very weak due to pair absorption, 
and below the detectability level for GLAST (unlike the case 
of Cyg X-l, see Zdziarski & Gierlihski 2004). 

Fig- shows a comparison of all the power spectra (in- 
cluding that of the ui state, see Section 15.11 for the entire 
PCA energy range, ~2-60 keV. All the x-state power spec- 
tra look rather similar, with no obvious dependence on the 
fitted strength of reflection. The power spectrum of VP 720 
is clearly the narrowest and shows the strongest QPO at the 
lowest frequency. The strong QPOs at frequencies of ~1.5- 
4 Hz are characteristic to the intermediate/very- high state 
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f [Hz] 

Figure 5. Comparison of the ~2-60 keV power spectra for both 
the x an d w states. 

of black-hole binaries (e.g., Zdziarski & Gierlhiski 2004). 
Also, the photon spectra (Figs. |2HHl an< i L/Le ~ 0.2 ap- 
pear characteristic to those of the intermediate/ very-high 
state (Zdziarski & Gierlinski 2004) This confirms this clas- 
sification of that state of GRS 1915+105 (Done et al. 2004) 
rather than as a hard/low state (Belloni et al. 2000). 

Fig. HJ shows the photon energy-dependent power spec- 
tra for three energy channels. We see that the power gen- 
erally increases with the photon energy. We also see that 
the width of the power spectra increases with the increasing 
photon energy. The main QPOs are surrounded by shoul- 
ders, showing weaker maxima at frequencies about a half 
and twice that of the QPO, with the low-frequency shoul- 
ders much more pronounced in the cases of VP 601 and 619. 



5.2 The oj state 

Fig. Q shows the unfolded spectra of the VP 813 uo state, 
corresponding to the high-level continuum and to the dips. 
The relative normalization of the OSSE spectrum is assumed 
to be unity (i.e., the same as of the PCA) for both broad- 
band spectra. This is due to having at our disposal only the 
integrated OSSE spectrum, with the relatively low statis- 
tics. The HEXTE data agree with the OSSE well. The two 
PCA/HEXTE spectra pivot at ~50 keV, and are thus con- 
sistent with the same level of the high-energy tail. 

The earlier version of the VP-813 data analyzed in Pa- 
per I required neither Compton reflection nor Fe Kq line. 
In contrast, we find that the present high-flux data do re- 
quire reflection, with the characteristic pattern of residuals 
in the fit without reflection with an edge at ~7 keV and a 
hump at > 10 keV. The probability of reflection appearing 
by chance is ~0.01, 0.004 for the compps and eqpair fits, 
respectively, using the F-test. Still, we cannot be sure it is 
not due either to some systematic effect in the data or to 
complexity of the X-ray absorber. Thus, we give in Table 2 
the results both allowing reflection and kept at null. On the 
other hand, the low-flux data do not require the presence 
of reflection, though the upper limits are relatively large, 



0,/2-k < 1.2, 1.9 in the COMPPS and eqpair models, respec- 
tively. 

Fig. Ufa) shows decomposition of the total intrinsic 
spectrum for the high-flux data as obtained by fitting the 
eqpair model with reflection (Table 2). The high-flux spec- 
trum and L/L-e ~ 0.5 are roughly similar to ultrasoft spec- 
tra of black-hole binaries (e.g., Zdziarski & Gierlinski 2004), 
although here strong Comptonization of the blackbody com- 
ponent, reflection, and relatively strong variability (Fig. 
see below) are required, in contrast to the standard ultrasoft 
state (Gierlinski & Done 2004). 

However, we have found that the decomposition of the 
prominent X-ray hump below ~20 keV into the unscattered 
disc blackbody and its thermal Comptonization is not com- 
pletely unique. We have considered a model with a sepa- 
rate (unscattered) disc blackbody component (diskbb) in 
addition to one with a hybrid plasma Comptonizing disc 
blackbody photons (eqpair) and allowing the two maximum 
blackbody temperatures of the direct disc emission and of 
the seed photons to the hybrid plasma to be different. This 
may correspond to a situation in which the plasma covers 
only an inner part of the disc while its outer part is seen 
directly. This model (hereafter abbreviated as DH) yields 
a number of relatively shallow local minima. In particular, 
we find a spectral solution with the maximum blackbody 
temperatures of the two components above of 1.48 and 2.43 
keV, respectively (the latter value requires the color correc- 
tion to be relatively large) at xt = 74/92, similar to the 
corresponding value of 76/94 in Table 2. The main param- 
eters of the hybrid plasma are £h/£ s = 0.39, £ n th/£h = 0.76, 
■n = 11.3, r inj = 2.4, Q./2-k = 0.7. The components of this 
model are shown in Fig. |HJb) . This model also required rel- 
atively strong relativistic smearing of Compton reflection 
(consistent with the model geometry), and thus we have 
fixed the inner disc radius at lOGM/c 2 . 

Although the fits with the hybrid model to the broad- 
band spectra are very good, we cannot rule out the origin of 
the high-energy tail from a separate process. We have thus 
fitted the model (hereafter abbreviated as TP) consisting of 
scattering of the disc blackbody emission by thermal elec- 
trons only (using eqpair with null reflection for the high- 
flux spectrum) and a separate power-law component. The 
obtained parameters for the thermal scattering of the disc 
blackbody are very similar to those given in Table 2, except 
that £h_/£ s is somewhat lower, 0.271q'q^ (since £h now does 
not include the nonthermal power). The fitted power law in- 
dex is F = 2.1±g |, fcT bb = 1.38l° °3 keV, n = 4.6±g f , and 
X l = 90/97. 

We note that the COMPPS fails to represent the high-flux 
spectrum because of the too high Thomson optical depth re- 
quired by the data. The iterative scattering method (Pouta- 
nen & Svensson 1996) used here utilizes only 50 scatter- 
ings, which number is much too low for the present case. 
The fitted value of r becomes then artificially very high. In 
our model in Table 2, we constrain r to 15, but this value 
is very strongly overestimated. Still, other parameters ob- 
tained using COMPPS (e.g., kT) appear correct, and they are 
in agreement with those obtained using eqpair. 

Fig. Qj] shows the comparison of the total PCA power 
spectrum of the entire VP-813 observation to those in the 
X state, and Fig. [!|]shows the energy-dependent PCA power 
spectra. We see the power spectra are distinctly different 
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Figure 6. The power spectra of the individual data sets corresponding to the energy ranges of ~2— 7 keV (red), 7—15 keV (green) and 
15-60 keV (blue). 



from those in the \ state. They have broad and flat (in 
power per logarithm of frequency) maxima at ~0.05-0.4 Hz 
and are cut off at both ends. 



6 FRACTIONAL VARIABILITY 

6.1 Data and multicomponent rms models 

Figs. llOl and llll show the results of calculating the rms using 
the two methods described in Section |3] for the binning of 
1/128 s. We see they yield very similar results. The small 
differences, with the PDS method yielding somewhat lower 
values of the rms at the lowest energies, are due to the fi- 
nite range of the frequencies integrated in that method. We 
note that the shown rms corresponds to most of the total 
variability power of the source, at least for / > IO -3 Hz, see 
Figs. El© |1 

We have modelled the rms dependencies by assum- 
ing contributions from the three fitted spectral components 
shown in Fig. [3] i.e., the unscattered blackbody disc emis- 
sion, Compton scattered component, and Compton reflec- 
tion. We assign each component a value of the fractional 
rms, ri(E) = ai(E)/Fi(E), which yields the total fractional 
variability of 



HE) 
F(E)' 



(3) 



in the case of the variability in the different components 
being completely correlated (coherent), and 



(4) 



in the case of uncorrelated (incoherent) variability. Here 
Fi(E) is the averege flux in the i-th spectral component, 
and F is the total flux. For notational simplicity, we omit 
hereafter the average sign for F, i.e., F(E) = (F(E)}, 
Fi = (Fi(E)). In general, r COI > r unc . 

We have tested that allowing the relative variability of 
the reflection, r rc fl, to be independent of the two other am- 
plitudes, the one for blackbody, r^b, and for scattering, r^b, 
only slightly improves the agreement between our models 
and the data. In fact, we expect the reflection component to 
respond to the incident flux from the Comptonizing plasma. 
Thus, we hereafter assume that the two components are fully 
correlated and r rc fl = r sc . 

We then first consider the case of both rt,b and r ac inde- 
pendent of energy. We have found this gives a fair descrip- 
tion of the rms observed during VP 601 and 720, as shown 
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Figure 7. Broad-band spectra in the u> state. The heavy (red) 
and and light (blue) curves give the spectra corresponding to the 
high and low flux states. The OSSE spectrum (green) is assumed 
to be at its own normalization at both states, which is consistent 
with the HEXTE spectra observed. The solid and dashed curves 
give the intrinsic and absorbed models, respectively. 

in Figs. HOt a). (c). The dashed and dotted curves show the 
cases of the correlated and uncorrelated variability, respec- 
tively. For VP 601, both give similar description of the data. 
However, those data were obtained in a mode allowing only 
very broad energy bins below ~15 keV, thus allowing only a 
rough test of the agreement. On the other hand, the data for 
VP 720 clearly prefer the correlated variability between the 
blackbody and scattered photons. (Note, though, that we 
have not performed here formal \ 2 fitting.) This is also the 
case for the data sets of VP 619 and 917, where the uncorre- 
lated variability predicts the rms spectrum more below the 
data than the correlated one. The cause of it is the dip in the 
rms profile around the intersection of the two components 
(each dominant at either low or high energies) due to par- 
tial cancellation of the two independent variability patterns. 
The coherence implied by the data (common in accreting 
black-hole sources, Vaughan & Nowak 1997) then can be 
understood by the feedback loop between seed blackbody 
photons and scattered photons. The Comptonizing plasma 
upscatters the seed photons, thus reacting to changes in their 
flux, but also the blackbody photons are partly due to the 
reprocessing of the emission of the hot plasma directed to- 
wards the disc. The presence of reprocessing is implied by 
the appearance of the strong Compton reflection in the data 
(see Table 2). 

The model with correlated, energy- independent, vari- 
ability in the two components gives a marginal description 
of the VP 619 (though the model curve lies below most of 
the data points), but it fails for the VP 917, see Figs. llOf b). 
(d). The disagreement occurs mostly in the region where the 
disc blackbody photons dominate the spectrum, see Fig. |3] 
Still, we find the flattening at the energy at which the un- 
scattered disc blackbody no longer contributes to the spec- 
trum roughly corresponds to the flattening in the observed 
rms dependence, which reassures us of the correctness of our 
general spectral decomposition. 
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Figure 8. Spectral components in the u-state spectra (solid 
curves). The dots show the unscattered blackbody, the short and 
long dashes show Compton scattering by thermal and nonther- 
mal electrons, respectively (the latter also includes the pair anni- 
hilation component), and the dot-dashed curves show Compton 
reflection. The panels (a) and (b) correspond to our main EQ- 
PAIR model with reflection (Table 2) and to the model (DH) with 
an additional disc blackbody component, respectively, see Section 
IP1 In (b), the dotted curve corresponds to the additional disc 
blackbody only, and the unscattered part of the emission of the 
Comptonizing plasma is not shown. 
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Figure 9. The power spectra of the lj state corresponding to the 
energy ranges of ~2-7 keV (red), 7-15 keV (green) and 15—60 
keV (blue). 



© 2005 RAS, MNRAS 000, 1-14 



GRS 1915+105: the distance, radiative processes and energy- dependent variability 9 



0.5 



0.4 



0.3 



0.2 



0.1 - 




10 
E [keV] 



50 



0.5 



0.4 



0.3 



0.2 



0.1 



incoh.: 
eoh.: r„, 
VP 720 



= 0.15, 0.24 
0.05, 0.24 



10 
E [keV] 



50 



0.5 



0.4 



100 



100 



incoh.: 
con.: r 



bb.max' ac' 

VP 619 



I >— 

„ = 0.15, 0.29 
= 0.11, 0.29 

= 0.55, 0.22, 0.06 




10 
E [keV] 



50 



100 



0.5 



0.4 



incoh.: r, 
coh 
r, 



t — | i i i i | 

r sc = 0.19, 0.33 
0.105, 0.33 



1 bb- 1 ac 

' h »™,. r„. r„ = 1.4, 0.295, 0.028 
VP 917 




10 
E [keV] 



50 



100 



Figure 10. The fractional rms variability as a function of energy for the x-state data. The heavy (blue) and light (red) crosses show the 
rms results obtained by the direct method and by integrating the PDS, respectively, for the binning of 1/128 s. The dotted and dashed 
curves show models based on the energy-independent contributions (as fitted by eqpair) from unscattered blackbody and scattered 
photons (including Compton reflection) in the correlated and uncorrelated cases, respectively, and with the parameters given on top of 
each panel. The solid curves in (a) and in (b), (d) assume radial variability of the rms in the scattered component and in the blackbody 
component, respectively, with the parameters given in the 3rd line of each panel. 



Thus, we have considered the rms of the disc blackbody 
component being energy-dependent. Physically, the emission 
in the disc-blackbody spectrum comes from a wide range of 
radii (Mitsuda et al. 1984), and the variability properties 
are expected to be dependent on the radius as perturbations 
propagate towards the black hole (Lyubarskii 1997; Kotov 
et al. 2001; Zycki 2003). The problem of rms from accretion 
flows has been studied by Zdziarski (2005, hereafter Z05), 
who has, in particular, obtained formulae for the rms of the 
disc blackbody assuming a power law radial dependence of 
the local rms. We have found that the model for rbb(-E) 
with variability of disc regions with the logarithmic length 
of Alni? = 0.1 each varying independently of each other, 
the value of fcTbb as in Table 2, the disc dissipation profile 
of Shakura & Sunyaev (1973) assuming the disc extending 
to the minimum stable orbit, and the radial dependence of 



the local variability index following the dissipation rate with 
the power law index of j3 = 4 [equations (11), (17) and (25) 
in Z05] gives good description to our data. (We stress again 
that we have not performed here a formal \ 2 fitting, and 
certainly there exist other models that can provide similarly 
good, or better, description.) The models, including the scat- 
tering component with the rms independent of energy corre- 
lated with the disc blackbody component, are shown by the 
solid curves in Fig. HUf b). (d). We have also found an addi- 
tional energy-independent rms component is required, which 
can be either correlated or uncorrelated with respect to the 
blackbody and scattering components. The cases shown in 
Fig. HUf b). (d) assume the E- independent rms component 
to be correlated with the other ones. 

We then see that the data for the VP 601 show the 
rms in the scattered component, > 15 keV, to continue to 
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increase. This may be due to a radial dependence for that 
component, similarly as for the blackbody component. We 
again used a model of Z05, assuming the local emission hard- 
ens from the photon index of 3.5 at 100 times the inner disc 
radius to 2.5 at the inner radius, the e-folding energy of 300 
keV, the fully correlated variability [equations (20-21) and 
(6) in Z05], and the dissipation profile as above except for 
(5 = 2. The resulting model is shown by the solid curve in 
Fig. 1101 Again, this is just one of many possible models, and 
the limited rms data at high energies prevent any better 
constraints. 

The rms data for the VP 813 (u state) are shown in 
Fig. 1111 We show here the rms spectrum corresponding to 
the entire data set, including the dips. However, the dips 
contribute relatively little to the total rms, and the rms for 
the high-flux state only looks almost the same. This is also 
related to the fact that the variability is the strongest in 
the high energy tail (Fig. 111! in spite of the spectra of the 
high and low-flux states (Fig. [7J differing mostly at low en- 
ergies. However, the transitions between those flux states 
occur only relatively infrequently, with most of the variabil- 
ity occuring within the high-flux state. 

In our main eqpair model either with or without reflec- 
tion (Table 2) , the flux in unscattered blackbody photons is 
comparable to the flux in scattered photons already at low 
energies, see Fig. |HJa) . Thus, assigning the scattered com- 
ponent even a very low energy-independent value already 
yields a substantial rms at low energies. Consequently, the 
upper limit on r sc is given by the rms at ~3 keV, and we re- 
quire rbb = 0. Such a model is shown by the dashed curve in 
Fig. [TT] and we see it fails completely at > 5 keV. Allowing 
the unscattered blackbody component alone to have a radial 
dependence, similarly to the cases of VP 619 and 917, does 
not improve the fit substantially. We also have considered 
the case in which the nonthermal electrons have a different 
variability amplitude than the thermal ones. However, this 
still does not improve the fit. 

On the other hand, we find that our alternative model 
with the additional disc blackbody component (DH), Fig. 
|5[b) in Section f5. 21 does provide a very good description of 
the rms profile. The dotted curve in Fig. 1111 shows a very 
good fit obtained by assuming the variability of the addi- 
tional disc blackbody and of the hybrid plasma being fully 
correlated, with the rms of the latter ^> that of the former. 

Since the hybrid plasma in this model covers an inner 
part of the disc, this model implies the local rms increasing 
towards the center, with the above two-component model 
being probably an approximation to a continuous rms ra- 
dial dependence. In order to test this possibility, we first use 
the values of Fi(E)/F(E) from the decomposition of the 
spectrum into thermally-Comptonized disc blackbody and 
a power law (model TP in Section fo,2t . We then assume the 
former varies as a disc blackbody with a power-law radial 
rms dependence (Z05). This is justified by the blackbody 
fit to the < 20 keV spectrum (which yields the maximum 
temperature of kTbb = 2.3 keV) approximating it to within 
10 per cent. We use the specific rms model of Fig. 9 in Z05, 
which uses the same assumptions (except for kTbb) as the 
models above for VP 619 and 917. The power law tail is as- 
sumed to vary independently of the main component. The 
increasing contribution of the power law at > 20 keV results 
then in a break in the rms dependence. This model yields a 
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Figure 11. The fractional rms variability as a function of en- 
ergy for the w state (VP 813). The data points and the dashed 
curve have the meaning analogous to that of Fig. [TO] The dot- 
ted curve shows a model corresponding to the fit including an 
additional disc blackbody component f Section 15. 21 . with the pa- 
rameters given in the second line. The solid curve shows a model 
based on the description of the spectrum by a disc blackbody with 
a local rms varying with radius and a power law, with the param- 
eters shown in the third line above. See Section 16.11 for details. 



very good description of the observed rms profile, as shown 
by the solid curve in Fig. 1111 We have found, however, that 
the data do not constrain the rms of the power law com- 
ponent, r p i, except for an upper limit of <0.3. Thus, as a 
minimum assumption, we set r p \ = 0. To account for the 
flattening of the rms spectrum at low energies, we also need 
an energy-independent overall variability, assumed here to 
be uncorrelated with respect to the blackbody-like compo- 
nent. 

6.2 Frequency dependence of the rms and energy 
spectra of the variable components 

In Section 16.11 we have used the rms integrated over the 
most (or entire) available frequency ranges, not taking into 
account the possible dependencies on a chosen frequency 
range. As seen in Figs. |S| and [5] the corresponding energy 
dependencies of the shape of the power spectra are rela- 
tively small for the \ state, and somewhat stronger for the 
u state. Detailed studies of those dependencies are beyond 
the scope of this paper. However, we do consider here the 
rms dependence for QPOs (present in the \ state). 

We first considered the VP-619 PCA data. We have 
fitted its energy-dependent power spectra by sums of 
Lorentzians (e.g., Nowak 2000; Belloni, Psaltis & van der 
Klis 2002), and calculated the QPO power from the am- 
plitude of its Lorentzian. Fig. 1121 shows the resulting QPO 
rms dependence together with the rms integrated over the 
1/512-64 Hz range (as in Fig.llOb'l. We see the shape of both 
dependencies are virtually identical at E > 5 keV, while the 
QPO rms decreases faster with the decreasing energy than 
the total rms at E < 5 keV. In the framework of our two- 
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Figure 12. The fractional rms variability as a function of energy 
for both the integrated PDS (the upper solid crosses) and the ~4 
Hz QPO (the lower crosses) of the VP 619. The dotted symbols 
show the rms of the QPO rescaled by a factor of 1.98. 



component, blackbody and Comptonization, model, this can 
be understood as the QPO variability following mostly that 
of the Comptonizing corona, and much less that of the black- 
body. 

A general feature of most of the rms dependencies 
shown in Figs. [MITTl is a flattening at E > 15 keV. This 
means that the spectrum of the variable component becomes 
identical to that of the time-average one. On the other hand, 
the variability is reduced at lower energies due to the pres- 
ence of a less varying component, presumably the disc black- 
body. 

This can be also seen from another representation of the 
energy dependent variability, namely the energy spectrum of 
a variable component, 



F var (£) = a(E) = r(E)F(E), 



(•») 



used, e.g., by Mendez et al. (1997) and Gilfanov, Revnivt- 
sev & Molkov (2003). Fig. I13f a) shows the variable spec- 
tra for the VP 720 corresponding to integration of the PDS 
in the ranges of 1/512-64 Hz, i.e., corresponding to most 
of the measured variability, and the 1.4-1.9 Hz, contain- 
ing the QPO. We see both variable spectra have relatively 
similar shapes, consistent with the result for the VP 619. 
They are also very similar to the shape of the average spec- 
trum at E > 15 keV, but fall below that at lower energies. 
This is again consistent with the picture in which the vari- 
able component, including the QPO, arises predominantly in 
the Comptonizing corona, whereas the blackbody emission 
is much less variable. The normalization of the 1/512-64 Hz 
spectrum at E > 10 keV, as shown by the lower solid curve, 
is 0.24 with respect to the average spectrum, i.e. equal to 
that used in the corresponding Fig. HUf c) for the rms of the 
Compton component. 

The dashed curve corresponds to the coherent super- 
position of the unscattered blackbody and Compton com- 
ponents with the relative normalizations as in Fig. HOf c). 
i.e., 0.05 and 0.24, respectively. We see that it approximates 




E [keV] 

Figure 13. The average, total variable, and QPO (from top to 
bottom) energy spectra for the (a) VP 720 and (b) VP 917. The 
solid curves give the best-fit models to the average spectra as 
well as those rescaled to the high-energy parts of the variable 
spectra. The dashed and dotted curve in (a) give models with a 
reduced unscattered blackbody component rescaled to the level of 
the variable and QPO spectra, respectively, see Section 16.21 The 
shown spectra have been obtained by fitting the EQPAIR model. 



well the variable spectrum, consistent with the good fit of 
the rms model in Fig. HOf c). We also see that the QPO spec- 
trum is very similar to that of the total variable component. 
It is well approximated by the dotted curve, which has the 
same shape as th dashed one, but is rescaled by a factor 2/3. 
Note that the exact shape of the deconvolved spectra in Fig. 
1131 is model-dependent to certain degree. 

Fig. I13f b) shows the analogous results for the VP 917. 
In this case, our fit to the data (Section does not allow 
a simple two-component modelling of the rms dependence 
(Fig. I10H ) , which may be due to the radial dependence of 
the rms of the blackbody (Section l6.H . Thus, we show here 
only the total model to the average spectrum renormalized 
to the levels of the variable components. We still see the 
same effects as before, i.e., the high-energy, now E>8 keV, 
parts of the variable spectra have the same shape as the 
average spectrum, and the QPO spectrum is approximately 
the same as the total variable one. 

We note that recently Rodriguez et al. (2004) have ar- 
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rived at a different conciusion. Nameiy, they ciaimed evi- 
dence for the presence of an additional component with weak 
variability at high energies. They based that conclusion on 
(i) an apparent decrease (fitted by an e-folded power law) of 
the QPO rms at E > 15 keV in a PCA observation of GRS 
1915+105 (obs. ID P80127-01-03-00) in the \ state and (ii) 
the corresponding PCA/HEXTE spectrum brighter at high 
energies than their fitted absorbed e-folded power law model 
(fig. 3 in Rodriguez et al. 2004). 

Given the potential importance of such a result for un- 
derstanding of the underlying physics of GRS 1915+105, we 
have repeated their rms analysis for the same observation. 
With the current PCA calibration, we have found only a 
flattening of the QPO rms(B) at high energies, very simi- 
lar to our results for the VP 619 and 917 above, with no 
evidence of any decline at highest energies. Also, that de- 
pendence is much better fitted by a broken power law (with 
the best-fit index ~ above the break) than by an e-folded 
power law (which improvement has the chance probability 
of 0.009 from the F-test). We also note that the rms decline 
at high energies was seen by Rodriguez et al. (2004, see their 
fig. 4) only in one out of 8 observations. In the remaining 
7 observations, a flattening at > 15 keV was observed, fully 
consistent with our results above. Furthermore, the e-folded 
power law used by them to fit the energy spectrum is clearly 
not a physical model. That spectrum is very similar to our 
X-state spectra (Fig.^, which are very well fitted by hybrid 
Comptonization, with no need for any additional component 
at high energies. 

Rodriguez et al. (2004) have also suggested that the pre- 
sumed additional high-energy component originated from 
the nonthermal synchrotron emission of a jet. We note here 
that even if there were evidence for such a component, this 
interpretation has severe problems. Namely, the high-energy 
part of the spectrum shown in fig. 3 of Rodriguez et al. 
(2004) is approximately of a power-law shape with F ~ 3. 
Such a steep component when extrapolated to low ener- 
gies would contain enormous luminosity (as pointed out by 
Zdziarski et al. 2003), as well it would greatly exceed the 
observed range of the IR flux, see, e.g., fig. 6 in Ueda et 
al. (2002). E.g., an extrapolation to ~ 10 15 Hz, where the 
nonthermal electrons would no longer cool efficiently (e.g., 
Heinz 2004), with T = 3 and at lower energies with F = 2.5 
(i.e., with the canonical AF due to radiative cooling) down 
to the turnover frequency at ~ 10 13 Hz yields the luminosity 
in that component as high as ~ 10 42 erg s _1 . A high-energy 
cutoff in a much harder nonthermal synchrotron spectrum 
fine-tuned to start at tens of keV would alleviate that prob- 
lem. However, the OSSE data show no evidence for such a 
cutoff in the \ state up to at least ~600 keV (Paper I) or 
more (Ueda et al. 2002). On the other hand, a nonthermal 
synchrotron cooling break in hard X-ray is highly unlikely 
(e.g., Heinz 2004). Summarizing, we can rule out the pres- 
ence of a nonthermal synchrotron component at hard X-rays 
in the \ state based on the spectral data alone. 

6.3 Alternative rms models 

An rms increasing with energy can be due to a pivoting 
process, in which the average spectrum is modulated pro- 
portionally to a varying power law fixed at the pivot energy, 
E p , for E > E p (Zdziarski et al. 2003, eq. [A6]). With E p w 2 



keV, this model with an addition of an _B-independent vari- 
ability can roughly (but significantly worse than the best 
models in Section 16. Ill describe the increase of the rms 
at E < 10 keV, but not the subsequent flattening. Also, it 
predicts the rms increasing with decreasing energy below 
E p . Pivoting occurs, e.g., when a thermally-Comptonizing 
plasma with an approximately constant power is irradiated 
by a variable flux of seed soft photons, which is characteristic 
to the hard state (Zdziarski et al. 2002). The pivot energy is 
then in the middle of the Comptonization spectrum, which is 
clearly not the case now. Thus, this model appears unlikely 
for the data presented here. 

A pattern of the rms increasing with energy occurs in 
the soft state of Cyg X-l on long time scales (Zdziarski et 
al. 2002). Those authors interpreted it as variability of the 
power supplied to the hybrid plasma at the power in the 
blackbody seed photons varying much less (Churazov, Gil- 
fanov & Revnivtsev 2001). This is similar to the models 
presented in Section [6. II More sophisticated models of this 
type could take into accout the changes of the spectral shape 
with the change of the plasma power. This could be mod- 
elled in the eqpair as variability of at the constant t B 
(Zdziarski et al. 2002). However, such detailed simulations 
are beyond the scope of this paper and will be presented 
elsewhere (Gierlihski & Zdziarski, in preparation). 

Yet another possibility is time evolution of flares. Zycki 
(2004, fig. 4), basing on the magnetic flare scenario of Pouta- 
nen & Fabian (1999), has found the rms increasing from 0.12 
at 1 keV to 0.21 at 30 keV. This increase is too shallow to 
account for most of our rms spectra, but these numbers de- 
pend on the model parameters used. On the other hand, it 
is unclear how this model can explain the flattening of the 
rms at high energies seen in our data. 



7 CONCLUSIONS 

We confirm the result of Paper I of a very good description 
of broad-band X-ray/soft 7-ray spectra of GRS 1915+105 
by hybrid Comptonization of disc blackbody emission, by 
using now all available RXTE/OSSE data in the \ an d w 
variability states. In the case of the u) state, we find that 
model applies to the spectra of both the high-flux steady 
emission and the dips. Analyzing the energy as well as power 
spectra, we identify the \ state with the intermediate/ very 
high state of black-hole binaries. In the u> state, we see strong 
and variable thermal Comptonization of the disc blackbody 
by a cool, Thomson-thick plasma, followed by a relatively 
weak high- energy tail. 

Our major new results concern the rms energy depen- 
dencies. We find the rms spectra can be well fitted by a 
coherent superposition of rms contributions from the fitted 
spectral components, namely the less-variable disc black- 
body and more variable Comptonization. Most of our data 
show a flattening of the rms at high energies, consistent with 
the domination of the Comptonization spectral component 
with the rms independent of energy. This is also the case for 
the rms energy dependence for the QPOs present in the x 
state, which have energy spectra similar to those of the total 
variable components. On the other hand, some of our data 
require the rms of the disc blackbody emission to increase 
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with energy, possibly due to the perturbation amplitude in 
the accretion flow increasing with the decreasing radius. 
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Table 2. Model parameters of the broad-band spectra. For each data set, the first model (COMPPS) fits the electron distribution, wheras the second (eqpair) fits the acceleration rate and solves 
for the electron distribution self-consistently. Two pairs of the results for the high-flux state of VP 813 correspond to fits with and without reflection, respecitvely. The quantities without stated 
errors are those derived by EQPAIR from the energy and pair balance (i.e., not model parameters), except for those denoted by 'f, which are parameters fixed during the fit. The ranges of Nn in 
both models are limited to > 1.8 X 10 22 cm -2 , and the range of r in COMPPS is limited to < 6 except for VP 813, where r is fixed at 15 (see Sections 4, 5.2). 
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